Phase separation of edge states in the integer quantum Hall 

regime 



L. Brey, J. J. Palacios, and C. Tejedor 

Departamento de Fisica de la Materia Condensada. Universidad Autdnoma de Madrid. 

CO ■ Cantoblanco, 28049, Madrid. Spain. 

ON 

on 



or 

< 

CO 
(N 



c3 



13 



Abstract 



Coulomb effects on the edge states of a two dimensional electron gas in 



J> . the presence of a high magnetic field are studied for different widths of the 

ON, 

boundaries. Schrodinger and Poisson equations are selfconsistently solved in 

O 

m . 

ON ■ appear for smooth interfaces in order to minimize the electrostatic energy re- 



tire integer Quantum Hall regime. Regions of flat bands at the Fermi level 



O , electron liquids. 

O 



lated to the existence of dipoles induced by the magnetic field. These plateaus 
determine the phase separation in stripes of compressible and incompressible 



A high magnetic field perpendicular to a two dimensional (2D) electron gas produces a 
spectrum with a continuous part in which states located close to the boundaries are involved. 
Those edge states are essential for the response of the system to external perturbations, par- 
ticularly in carrying a net electrical current when a bias is applied to the sample |IJ . For an 
empty system (or when electron-electron interaction is neglected), edge states are completely 
determined by the confining potential. However, when a density of electrons populates the 
system, Coulomb interaction changes qualitatively the band structure. Such a problem has 
been only addressed in the case of a very smooth external potential either by a Thomas- Fermi 
approach ||[J or by solving the Poisson equation in the classical limit Spin splitting 
effects have also been studied in a variational scheme [||. In some cases, like in the case of 
the Fractional Quantum Hall regime, in which many-body effects are crucial, only general 



arguments have been given to discuss the importance and behavior of edge states |6|-|T0||. In 



all those works, the smoothness of the confining potential is a necessity due to the strong 
simplifications involved. The actual situations in experiments cover the whole range from 
abrupt (etching techniques) to smooth (gate techniques) boundaries (T|]. To understand the 
properties of the edge states for any width of the interface region a quantum analysis in- 
cluding electrostatic effects is needed. This implies to solve selfconsistently the Schrodinger 
and Poisson equations for electrons confined by a potential defined in a 2D system in the 
presence of the external magnetic field. 

We are interested in studying the properties of electrons confined in the z-direction by 
semiconductor interfaces or quantum wells affected by a magnetic field B = Bu z . The system 
has no constrictions in the y-direction while in the x-direction there are some boundaries 
with typical widths significantly larger than the extent of wave functions in the z-direction. 
Therefore, we neglect the effects of the width in z of the electron gas and just consider a 
strictly 2D (xy) system with boundaries in the x-direction. Physically, such boundaries are 
produced by gate potentials which deplete the 2D gas in some regions confining the electrons 
in the rest of the xy plane. We have a total depletion under the gates by putting infinite 
barriers in x = and x = W, while the smoother interface region is the sum of two terms: 
the electrostatic potential produced by the electronic charge, and a confining potential in the 
x-direction which we take as the one created by a fictitious distribution of positive charge 
with the trapezoidal shape 



< \ n + 



x < x < L 

L L<x<W-L (1) 

W-x W - L < x <W 



In this way, the electrons with a 2D density n + are smoothly confined in a channel of width 
W . With this model, the Landau gauge A = (0, Bx, 0) is the adequate to study the problem. 
The wave functions have the form 



2 



e iky 



* n ,k(x,y) = -f=<l>° k (x) (2) 



where the wave function 0^ k ( x ) * s an eigenstate, with eigenvalue e a n k , of the one-dimensional 
Schrodinger equation 



h 2 d 2 m*u 2 (x - x ) 2 , ag*p B B 



+ C ~\ — + eV(x) + 



= (3) 



2m* dx 2 2 v ' 2 

where u; c = eB/m* is the cyclotron frequency, xq = k£ 2 is the semiclassical center of the 



orbit in terms of the magnetic length £ = y (K/eB), a is the spin of the electron, g* is the 
effective ^-factor and fiB the Bohr magneton. The selfconsistent electrostatic potential is 
given by 

V(x) = — f dx'[p + {x') - p{x')} In \x - x'\ (4) 



e being the dielectric constant of the medium and the electronic selfconsistent charge density 
being given by 

rEp 

p(x) = J2 del k \r ntk {x)\ 2 (5) 

cr,n J 

Ep is the Fermi level which, for given total charge density, must be computed within the 
iterative procedure to solve the Schrodinger equation. We define the total charge by means 
of the filling factor v for the 2D bulk (i.e. the centre region of the wide channel). The iter- 
ations are started by taking the electronic charge locally equal to the positive background, 
i.e. with a potential V(x) which is flat in between the two infinite barriers. We neglect the 
exchange interaction between electrons, therefore we restrict our calculation to even filling 



factors where exchange effects are less important | IT| . Electrostatic effects are the same for 
any v corresponding to an incompressible liquid in bulk, i.e. with a large enough gap at the 
Fermi level. Therefore, it is possible to extend our conclusions to odd filling factors. 



Since it is very interesting to discuss what happens to the different branches of the disper- 
sion relation, we present results for the case v — 6 and B = IT. We use typical parameters 
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for GaAs so that g* is so small that there is no difference between spin-up and spin-down 
bands. In order to analyze the importance of the boundary smoothness figures 1 to 3 give 
the dispersion relations e n ^ = e n (xo), potential profiles V(x) and charge densities p(x) and 
p+(x) — p(x) (multiplied by 50 ) for the cases L = 10£, L = 30£ and L = 50£ respectively 
All the figures give results computed with W = 2L + 10£ in which the two interfaces are 
decoupled from each other as we have checked out by increasing the separation (changing 
the factor 10 by higher numbers) and obtaining the same results. Due to the symmetry with 
respect to the centre of the channel we only present results for x < W/2. From figure 1, one 
observes that an interface of ten magnetic lengths can be considered rather abrupt because 
no plateaus at the edge appear in V(x) or e n (x ) as suggested for smooth interfaces J4|. Only 
some inflections appear for L = 10£ that become flat regions for L = 30£ and still broader 
for the smoother case of L = 50£. For L = 30£ and L = 50£ there are two types of plateaus 
near the Fermi level which are intrinsically different from each other. Those spatially close 
to the edges are not completely occupied because they lie in energy exactly at the Fermi level 
while the plateau in the bulk is clearly below Ep so that its states are fully occupied. Since 
we work in a single particle model, our results for the flat regions at the edge present very 
small oscillations instead of being completely flat so that they accommodate both electrons 
and holes around the Fermi level. In our calculations, the amplitude of such oscillations 
tends to zero, and the plateaus can be considered absolutely flat for any physical purpose, 
in particular for considering these regions as locally compressible liquids corresponding to 
noninteger filling factors. 

Let us try to understand the physical causes for the shape of the dispersion relation by 
starting with the flat regions close to the edge. Among these plateaus, there are regions 
of incompressible liquid corresponding to integer local filling factors. When L increases, 
the compressible regions extend while the incompressible ones remain practically the same. 
This observation helps to understand the physical origin of the spatial separation of the 
compressible and incompressible phases. In a region of integer filling factor, the electronic 
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charge density is practically constant and the positive background varies linearly. This 
implies an electrostatic dipole that increases with the width of the incompressible region. 
Since that dipole means an increase of the energy, the system reacts to minimize the energy 
by decreasing the width of the incompressible region as much as possible i.e. by forming 
compressible regions. In this new phase, the electronic charge does not need to be constant 
any more. Therefore it can have the same shape of the positive background (as shown in 
the figures) eluding any increase of the electrostatic energy. It is important to stress that 
this argument is independent of the form of a background which varies with the position so 
that our results should be completely general for any interface. In the bulk of the sample, 
the fictitious positive background is constant, the Fermi level lies above the Landau state 
giving the well known incompressible liquid for integer filling factor. The width of the in- 
compressible phases at the interface can not decrease up to zero because there is a lower 
limit imposed by the extent of the wave functions. The minimization of the electrostatic 
energy, which depends on the strength e 2 jel of the interaction, dominates on any variation 
of kinetic energy, which depends on the cyclotron energy Hu c . 



Disorder induced localization only affects to the bulk electronic structure but not to edge 
states so that our results are valid for actual samples 1^]. The experiments for integer v 
made on samples with boundaries created by etching involve abrupt interfaces and the elec- 
tron liquid must be incompressible. On the contrary, for samples with boundaries defined by 
gate potentials the interfaces are typically broader than ten times the magnetic length and 
phase separation should occur. The current will be carried in the compressible regions where 
the zero velocity of each state compensates with the infinite density of states. So, transport 
does not seem to be a good way to detect the electrostatic effects here studied. Spectroscopic 
experiments are better candidates to analyze compressible regions because the high density 
of states must produce strong Fermi edge singularities in absorption and emission of light 
as well as significant alterations in the spectrum of edge magnetoplasmons. 



5 



This work has been supported in part by the Comision Interministerial de Ciencia y 
Tecnologia of Spain under contract MAT 91 0210 and by the Commission of the European 
Communities under contract SSC - CT 90 - 0020. 



6 



REFERENCES 

For a review see e.g. C. W. J. Beenakker and H. van Houten, in "Solid State Physics" 
ed. by H Ehrenreich and D. Turnbull (Academic, New York, 1991), Vol. 44. p. 1 

P. L. McEuen et al., Phys. Rev. B 45, 11419 (1992). 

I. K: Marmorkos and C. W. J. Beenakker, Phys. Rev. B 46, nov. 92. 

D. B. Chklovskii, B. I. Shklovskii and L. I. Glazman, Phys. Rev. B 46, 4026 (1992). 

J. Dempsey, B. Y. Gelfand and B. I. Halperin, preprint. 

C. W. J. Beenakker, Phys. Rev. Lett. 64, 216 (1990). 

A. H. MacDonald, Phys. Rev. Lett. 64, 220 (1990). 

X. G. Wen, Phys. Rev. Lett. 64, 2206 (1990). 

A. M. Chang, Solid St. Commun. 74, 871 (1990). 

H. A. Fertig and S. DasSarma, Surf. Sc. 263, 60 (1992). 
C. Kallin and B. I. Halperin, Phys. Rev. B 30, 5655 (1984). 

F. J. Rodriguez and C. Tejedor, to be published in Phys. Rev. B 47, 15 january (1993) 
and references therein. 

B. I. Halperin, Phys. Rev. B 25, 2185 (1982). 



FIGURES 

FIG. 1. Dispersion relations e n (xo) (broad continuous lines), potential profiles V(x) (interme- 
diate continuous line) and charge densities p(x) (narrow continuous line) and 50[p+(x) — p{x)\ 
(dashed line) as a function of the position for a boundary with L = 10£. The Fermi level (dotted 
line) selfconsistently computed is also shown. 

FIG. 2. Dispersion relations e n (xo) (broad continuous lines), potential profiles V(x) (interme- 
diate continuous line) and charge densities p{x) (narrow continuous line) and 50[p+(x) — p(x)] 
(dashed line) as a function of the position for a boundary with L = 30£. The Fermi level (dotted 
line) selfconsistently computed is also shown. 

FIG. 3. Dispersion relations e n (xo) (broad continuous lines), potential profiles V(x) (interme- 
diate continuous line) and charge densities p(x) (narrow continuous line) and 50[p+(x) — p(x)} 
(dashed line) as a function of the position for a boundary with L = 50£. The Fermi level (dotted 
line) selfconsistently computed is also shown. 
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